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Abstract 

Methods of Pade approximation are used to analyse a multivariate 
Markov transform which has been recently introduced by the authors. 
The first main result is a characterization of the rationality of the 
Markov transform via Hankel determinants. The second main result 
is a cubature formula for a special class of measures. 
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1 Introduction 

Let (T be a non-negative finite measure on a subinterval [a, b] of the real line 
R. Then the numbers J x^da [x) are called the moments of the measure a. 
The Markov transform of a is defined for C G C \ [a, b] by the formula 



^(0:= / ^da{x). (1) 
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In the theory of moments Pade approximation of the Markov transform a (C) 
is an important tool, see pQ, [S], [B] or [TS] and section 6. Here Pade approxi- 
mation is performed at the point oo, so we consider the asymptotic expansion 



^(C) = E / x'daix)—^ioT \C\>R. (2) 

Let now /x be a signed measure fj, on the euchdean space M'' with support 
in the closed ball Br := e M'^ : < R} where |x| := + ... + is 
the euclidean distance for x = {xi,...,Xd) G M"^. In [2] we introduced a 
multivariate Markov transform for the measure n by the formula 

'^^^'^) = / ^:^^^/^(^) Kl > e §^-^ (3) 

Here S"'^^ := |x G M*^ : |x| = l} is the unit sphere, and C is a complex number 
with ICI > R. In the denominator, the expression r ((^^ — x) is the analytic 
continuation of the function p i — > \p9 — x\ defined for p G M with p > R, see 
Section 3 for details. The motivation for this definition stems from the work 
of N. Aronszajn about polyharmonic functions and the work of L.K. Hua 
about harmonic analysis on Lie groups, see [2], [11] or [T3]. Following the 
analogy with the one- dimensional case, we consider the asymptotic expansion 
of the multivariate Markov transform. From the growth behaviour at infinity 
of the kernel C'^^^ jr {C^Q — x^ it is easily seen that the asymptotic expansion 
is of the form 

oo 

/^(C^) = E/'(^)7m (4) 

for \(\ > R and 9 G S'^"^ where // : E>'^~^ — )■ C are continuous func- 
tions. The aim of this paper is to show that methods from Pade approx- 
imation can be successfully used for an analysis of the multivariate Markov 
transform. Roughly speaking, we shall perform in (jl]) the classical uni- 
variate Pade approximation for each fixed 6 G S^^^ obtaining a Pade pair 

Let us describe the results in the paper: In section 2 we shall first review 
the basic notions from Pade approximation which are needed in the paper. 
In section 3 the asymptotic expansion defined in will be investigated. It 
turns out that each coefficient function fi in (jl]) is a finite sum of spherical 
harmonics of degree < I, and each fi is the restriction of a homogeneous 
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polynomial Fi (x) of degree / to the unit sphere. The Hankel determinant 
of the multivariate Markov transform /i (or a measure /i) is defined by the 
expression 

//oW fi{0) ■■■ fn-i{e) \ 



Hn (/i, 6) := det 



(5) 



V fn-l (0) fn (0) ■■■ /2n-2 (0) J 

In section 4 we show that the Hankel determinant iJ„ {fi, 9) of a measure is 
the restriction of a homogeneous polynomial of degree n (n — 1) to the unit 
sphere. In section 5 we shall prove a Kronecker type theorem: the Hankel 
determinants Hn (/i, 9) are zero for all large n if and only if the function 
C, I — )■ Jl{(,9) is rational for each 9 G S*^"^. Moreover, this is equivalent to 
the rationality of the multivariate Markov transform Ji. 

A measure /i is called Hankel positive if the Hankel determinants Hn {fi, 9) 
are strictly positive for all natural numbers n and for all 9 G S'^~^. In section 6 
we prove that for each Hankel positive measure fi there exists a non-negative 
measure fin which is equal to fi for all polynomials of degree < 2n — 1 and 
which has support contained in an algebraic variety. Further we character- 
ize Hankel positivity by an extension property of the multivariate Markov 
transform. 

Finally we need some notations from harmonic analysis. A function Y : 
E)'^^^ — )■ C is called a spherical harmonic of degree G No if there exists 
a homogeneous harmonic polynomial P (x) of degree k (in general, with 
complex coefficients) such that P (9) = Y (9) for all 9 G S"^"^. Throughout 
the paper we assume that Y^^rn (x), m = 1, a^, is a basis of the set of all 
harmonic homogeneous polynomials of degree k which are orthonormal with 
respect to scalar product 



{f,9)s^-. ■■= / f{9)g{9)d9. 

Here denotes the dimension of the space of all harmonic homogeneous 
polynomials of degree k. By cua we denote the surface area of S'^~^. 



2 Basic facts from Pade Approximation 

At first let us recall some basic facts from Pade Approximation (we refer to 
[IB] for proofs): let / be a holomorphic function for ^ G C, |C| > -R, of the 
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form 



oo ^ 



1=0 



Let n be a natural number. Then there exists a polynomial P„ 7^ of degree 
< n such that 



(6) 



l=n 



where Qn is the polynomial part of the series P„ (z) / (C) ; it is easy to see 
that Qn has degree < n — 1. A pair (P„, Q^) is called an n-th Fade pair if P„ 
and Q„ are polynomials, P„ 7^ 0, degP„ < n and degQ„ <n — l, and they 
satisfy ([6]). An index n is called normal if for any n-th Pade pair (P„, Qn) the 
polynomial C, ^ Pn (C) has degree exactly n. Proposition 3.2 in [18] shows 
that n is normal if and only if the Hankel determinant 



//o 

/i 



fi 

f2 



\ /n— 1 /n 

is not zero. If n is normal then the polynomial 

/ /o /i 

Pn (C) := det 



fn 

f2n-2 ) 
fn 



(7) 



fn-l /2n-l 

VI c ■ ■ ■ 

has exact degree n and (P„, Qn) is an n-th Pade pair where Qn is the poly- 
nomial part of Pn {z) f (z) . For arbitrary n, the rational function 

Qn (0 



(C): = 



is called the n-th diagonal Pade approximant of /. 



3 Asymptotic expansion of the multivariate 
Markov transform 

Following the exposition in [21 Section 2.2] we show that the multivariate 
Markov transform is well-defined. Let us set r (x) := . For p > and 
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6 G and x = {xi, x^) we have {p6 — x) = — 2p {6, x) + \x^ where 
{9, x) is the usual inner product in R'^. We replace p by a complex number C, 
and obtain 

r2 [ce -x) = c^-2C {e, x) + = (c - x)f + \x\^ - I (0, x) 1^ . 

Note that \x\^ - \ {B,x)\^ > for each 6 G S'^"^ If we define 



a {e, x) := x) + i^\xf - \{e,x)f 

then (C^ - x) = (C-a(^,x)) (^C - • Since |a(^,x)|^ = |x|^ it fol- 

lows that ((^^^ — x) 7^ for all \(\ > \x\ . Next we see that the function g, 
defined by 



^^^y_ rHce-x) _ / a{e,x) \ f a{e,x) 



e \ c J \ c 

for \(\ > \x\ , has the property that g {() ^ (— oo,0]: since 



a{e,x) 

c 



< 1 and 



a(e,x) 

c 



< 1 it follows that 1 — 1 — ^(^2£) qj-q the right half plane, 



i.e. that their real parts are strictly positive, and therefore g {() ^ {—oo, 0] . 
Using the square root function a/^ defined on C \ (— oo, 0] one can define for 
\(\ > \x\ the analytic function 



c 



'r2 (C^-x) 



It follows from these facts that the multivariate Markov transform p (C, ^) is 
well-defined. 

Further we will make use of a real version of the multivariate Markov 
transform which we define by (note that we use d instead of ci — 1 as exponent 
in the nominator) 

Arcai (y) ■■= / — ^rfp (x) for yeR'' with \y\ > R. (8) 

jRd r[y — x) 

The real Markov transform Preai (z/) is related to Ji{(,6) in the following 
way: using results about harmonicity hulls and Lie norms (see P- 64]) one 
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may show that the function y i — V /ircai {v) has an holomorphic extension to a 
natural set Cr in the complex space C^, and the extension will be denoted by 
Aureal {z) for complex z G Cr. The set Cr is the set of all z = (zi, Zd) G 
such that 



L^{z):=^j\z\'-^J\zt-\q{z)\'>R (9) 

where we have defined \z\^ = \zif + ... + \zdf and q (z) = zf + ... + The set 
Cr is connected and open, and it contains all points ( ■ 9 with ( ^ C, \(\ > R 
and 6 G E>'^~^. The Markov transforms /I (C, 6*) and /ireai (-2) are related by the 
simple formula 

/"real (C^) = C/" (C, ^) for all C G C, |C| > i? and ^ G (10) 

Next we want to describe the asymptotic expansion of the multivariate 
Markov transform. Using the Gaufi decomposition of a polynomial (see The- 
orem 5.5 in [3], [22], or [13]) it is easy to see that the system 

kl^" yfc,m (x) , s, G No, m = 1, afc 

is a basis of the set of all polynomials. The numbers 



Cs,k,m-= I \x\'^'Yk^^{x)diJ,{x), s. A; G No,m = l,...,afc (11) 



are sometimes called the distributed moments, see [I2] . For a treatment and 
formulation of the multivariate moment problem we refer to |9] and [2]. From 
we cite 



Theorem 1 Let fi be a signed measure on M'^ with support in the closed ball 
Br. Then for all \C\ > R and for all 9 G E>'^^^ the following relation holds 

00 00 <^k -tr / n\ f 

/^(c,^) = EEE7S^ / \A''y^)d^^{^)■ (12) 

„_n n 1 



=0 fc=0 m-- 



For // defined in (jl]), a rearrangement of the series f|T2l) in powers C'"*"^ 
yields the relation 



[I] ai^2t 



fl (^) = J2Y1 (^t,l-2t,n.Yi^2t,m (9) , (13) 
(=0 m=l 

where [x] denotes the largest integer n such that n <x. 
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Proposition 2 For each I G Nq the coefficient function fi in is a finite 
sum of spherical harmonics of degree < I. Moreover, there exists a homoge- 
neous polynomial Fi (x) of degree I such that 

Fi {(6) = Cfi (6) for all 6 e S'^'^ and C G C. 

Proof. Formula f[T^ shows that // is a sum of sphercial harmonics of 
degree < I. Define a homogeneous polynomial Fi of degree I by 

[5] Oll-2t 

Fi (x) := XI "^M-si,™ Yi_2t,m (x) . (14) 

t=0 m=l 

By inserting x = p6' in (HM for positive p we obtain Fi {p6) = fi {6) . Since 
p I — y Fi {p6) is holomorphic we may replace p by a complex number (. The 
proof is finished. ■ 

The coefficient function /; can also be described by Legendre polynomials 
Pk (t) of degree k and dimension d, for definition see [17]. Clearly ( !T3|) and 
(ITT]) implies that 

„ [1/2] ai_2t 
fl(^)= / 5Zk|'*5^yi-2i,m(x)-yz_2i,r„(^)rf/x(a;). 

J^" t=0 m=l 

The addition theorem for spherical harmonics (see [T7]) says that 
Yk,^ (x) ■ n,™ {e) = \x\'' akPk f/^, , 

m=l \ \ I I / / 

SO one obtains the alternative description 

fi (e) = J2 «'-2t J^^ \4 Pi-2t (^^^, OJj dp (x) . 

We conclude this section with some examples and results illustrating the 
definitions. 

Example 3 Let a he a finite non-negative measure on an interval [a, h] with 
a > and consider the measure p = a^d6, i.e. for every continuous function 
f holds 

f{x)dp:= [ [ f {r9)da{r)d9. 



a JS' 



Then the distributed moments Cs^k,m cife zero for all k > since Yk,m (^) 
orthogonal to the constant function with respect to the measure dO. Hence 
[W\) shows that 



A measure /i on M*^ is called rotation invariant if /i (T^^ (5)) = /i (B) for 
all Borel sets B and for all orthogonal linear maps T : R'^ — J-M'^. The follow- 
ing result shows that a rotation invariant measure has a Markov transform 
^{(,,6) which does not depend on G Ef^^^ . Since the result is not needed 
later we omit the proof. 

Theorem 4 Let ^ be a measure on with support in Br. Then fi (C, 6) is 
independent of 9 if and only if fi is rotation invariant. In that case the mul- 
tivariate Markov transform possesses an analytic continuation to the upper 
half plane, namely 



for all ImC>0 andO e 

4 Multivariate Pade approximation and Han- 
kel determinants 

We start with the following observation: 

Proposition 5 Let Hn (/i, 9) be the Hankel determinant defined in Then 
there exists a homogeneous polynomial Hn {x) of degree n{n — 1) such 




From this we conclude that for all I G No and G S' 





(/i, 9) for all 9eS 



id-l 
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Proof. By Proposition |2] there exists a homogeneous polynomial Fi (x) 
of degree I such that Fi {(9) = (^fi {6) . Let us define 



Hn (x) := det 



Fi (x) F2 (x) 
\ (x) Fn (x) 



F„-i(x) \ 

F2n-2{X) J 



Now we replace x by ^6' and we apply the Leibniz formula for determinants to 



the matrix A 
Then 



[a. 



defined by a^j = Fi+j {(6) for i,j = 0, ...,n- 1. 



sign (a) Fo+<^(o) (C^*) ...Fn-i+a{n-i) {(d) ■ 

a permutation 

Note that 

+ a (0) + 1 + a {1) + ... + {n - 1) + a {n - 1) = n {n - 1) . 

It is obvious that Hn (x) is a homogeneous polynomial of degree n{n — 1) . 
We can factor out ("("-i) and we see that {(0) = C^''-^^Hn (/i, 9) . ■ 

In the following it is convenient to introduce the following notation: for 
a natural number n define a polynomial P„ (C, 0) of a univariate variable C 
of degree < n by 



PniC^e) :=det 



/ /o (^) /i (^) 



/2n-l (^) 



(15) 



We shall also write 

Pn (C, 0) = Po (e) +p,ie)c + ... + Pn (e) C- 

We define Qn (C? ^) as the polynomial part of Pn {(, 0) /i (C, 0) , so 

Qn (C, 0) = PnfoC"' + (Pn-l/o + Pnfl) C"' + - + (Pl/o + P2/1 + •.• + Pn/, 

From the results in section 2 the following is clear: 



(16) 



n-l) 

(17) 
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Theorem 6 // Hn (/i, 0)^0 then (P„ {(, 9) , Qn (C, d)) is an n-th Fade pair 
for the function 

oo ^ 

C^Ji{Ce) = J2fd0)-^, 

1=0 ^ 

for \(\ > R where 9 G S^^^ acts as a parameter. 

In Example |9] below we shall show that Q i — > Pn (C? ^) i^^-Y be the zero 
polynomial for certain 9 G S'^^^, so {Pn (Cj ^) ? (Cj ^)) is not always an n-th 
Pade pair. 

The advantage of working with P„ ((^, 6') is seen from the following result: 

Theorem 7 Let P„ (C, 9) and Qn (C, ^) be defined in ( f73]) and ([I^. T/ien 
i/iere exists a polynomial An of degree <n^ + n and a polynomial Bn of degree 
< + n — 2 such that 



C Pn (C, 0) = An m and C Qn (C, 0) = (Bn {(9) . 



(18) 



Proof. By Proposition [2] there exists for each / G No a homogeneous 
polynomial Fi of degree / such that (^) = Fi {(9) . Let us multiply each 
j-th column in (fT5|) . j = 0, n, with Let us define dn to be the sum 

of n - 1 + j for j = 0, n. It follows that C'^'^Pn {(, 9) is equal to 



det 



I2n-1 
'•271— 1 /-n 



C-'fn-i{e) ■■• 

Since Fi {(9) = C} fi {9) we obtain that C^^Pn (C, d) is equal to 

/ c~'Fom c~'F^ico) ■■■ c-'Fnm\ 

"^^^ Fn-i{C9) ■■■ ■■■ F^n-im 



From the j-th row factor out ^ ^ for j = 0, n — 1 and from the last one 

-n-l rpV,OTi rl _ /'^ _ 1 ^ _ Y^'^-l 

^j = 



C-\ Then rf„ 

C'P„(^,^)=det 



n — 1) — X]j=o ^ i^ equal to n^. Hence we have proved that 
/ FoiCO) F,m ■■■ FniCO) \ 



Fn-im ■■■ 

\ 1 



F2n-l m 
^2n 
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It follows that C Pn {(,0) = An {(6) where A„ (x) is defined as 





/ 


Foix) 


Fi{x) ■ 


■■ Fnix) 


\ 


An (x) := det 
















Fn-l {X) 




■ ■ F2n~l {x) 






v 


1 


1 |2 

F 


1 |2n 

\x\ 


1 



(19) 



This formula shows that the degree of A^ (x) is lower or equal than -n? + n. 
Let us discuss the polynomial part Qn (C, ^) • Let us write P„ 6*) = 
(^) +Pi (6') ^ + ... +p„ (6') By formula (IT3i) it is clear that (6') can be 
defined by the determinant of the matrix in f|T5|) where we have deleted the 
j-column and the last row. An analysis analog to the above shows that there 
exists a polynomial Rj (x) such that Rj {Qd) = ~^pj {6) . Now formula ([1] 
shows that 



n—l n—l—k 



k=0 



1=0 



Since^C"'CVi iO)Pk+i+i (0) = C'+^Rk+i+i {(9) Fi {(0) one can conclude that 
^C" Qn (C) 9) is a polynomial. ■ 

We want to relate the Pade approximation in Theorem [6] to Fade ap- 
proximation in the context of polynomials in several real variables. Let Fi 
be the homogeneous polynomial of degree / defined in Proposition |2l The 
asymptotic expansion 



oo ^ oo ^ 



1=0 



1=0 



and the identity yUrcai (C^) = C/^lC;^)? see ffTOl) . yield the asymptotic expan- 
sion of the real Markov transform /ij-eai iu) ; namely 



oo ^ 
/^real(2/) = y]^K^/) 

By formula Theorem [6] and [7] we can find polynomials An {y) and Bn {y) 
such that 

oo ^ 

An m /2 (c, 9) - (Bn m = E ^ 



l=n 
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for all 6 E E>'^ ^ such that the index n is normal, i.e. Hn {fJ^, 9) 7^ 0. We 
multiply this equation by C and write 

An m Cfi (c, e) - m = J2 (c^) t^- 

l=n ^ 

Further for the polynomial h{y) = \yf we have h{(6) = C^, so the last 
equation implies for real y = pO with \y\ > R and Hn (/U, 6^) 7^ 

An (y) /ireal (v) " Ivf Bn (y) = \y\'' V Fi (y) 

i=n \y\ 

Here An (y) and Bn (y) are subject to the conditions expressed in (ITS]) , and 
it seems to be rather technical to convert this in direct conditions for An, Bn- 

We refer to [2j and |8j for multivariate Fade approximation based on 
polynomials in several variables. 

Now we want to give an example of a measure /i such that the polynomial 
( I — > Pn (C, 6*) (defined in ( ITSl) ) is the zero polynomial. We recall at first the 
following result from [H]. 

Proposition 8 Let a be a measure on M with compact support, 6q be the 
Dirac measure on M at the point and let p = a ® 5q be the product measure. 
Then the multivariate Markov transform /i is given by 

- — / ^ H\ 1 ^--^ /" / , / X sin (/ + 1) t 1 . . 

a ® 5o (C, e^*) = - E y ^ (x) ^ (20) 

where U2 is the area measure o/S^. 

The last proposition has been used to show that there exists a measure 
/i with a support contained in an algebraic set such that (fi {(, 0) is not a 
rational function. 

Example 9 Let a be the Lebesgue measure on [0,1], so j^x^da{x) = l/(/ + 
1) and let n = a ® 5{) as in Proposition^^ Then // (1) = 1 for all I eNq, and 
this implies that Pn (C, 1) = for all n> 2. 



12 



5 Rationality of the multivariate Markov trans- 
form 



Recall that a function / : M'^ — )■ C is rational if there exist polynomials p {x) 
and q (x) with f (x) = p (x) /q (x) for all x with q (x) ^ 0. 

A theorem of Kronecker (Theorem 3.1 in [IB]) says that a necessary and 
sufficient condition for a series / {() of a single variable ( to be the Laurent 
expansion of a rational function is that the Hankel determinants Hm (/) are 
zero for all sufficiently large m. 

We have now the following analogue for the multivariate Markov trans- 
form: 

Theorem 10 Let fi be a measure with support in B^. Then the following 
statements are equivalent: 

a) For each 9 G 'B'^^^ the function ( ^ Ji{C, 9) is rational. 

b) There exists n G N such that {fi, 6) = for allm > n and 9 G S*^"^. 

c) There exists ri G N such that for each 9 G the function ( ^ Ji{C, 9) 
is rational of degree < n. 

d) There exist polynomials P {x) and Q (x) such that for all 9 G E>'^^^ and 
for all Id > R 

Pm^^andfi{U) = Cj^y (21) 

Proof. Assume a) and let d (9) be the degree of the rational function 
( i-> Ji{(,9) for 9 G S'^^^ (recall that the degree of a rational function / = 
p/q with relatively prime polynomials is defined as max {degp, degg}, see 
[m p. 38]). By Kronecker's theorem (cf. [HI p. 46]) it follows that the 
Pade approximant 7r„(C,6') is equal to 12 {(,9) for all indices n> d{9) and 
Hn (/U, ^) = for all n > d {9) . It follows that S"^"^ is the union of the 
following sets 

A„ := n^=„{0G§^-i:i7^(/x,^) = O} 

for n G N. Moreover is closed by continuity of 6^ Hn{fi,9). By 
Baire's category theorem there exists an index n such that An contains 
an interior point. Hence there exists 9q G S"^"^ and a neighborhood U 
of 9q such that Hm (/i, 9) = for all 9 ^ U and for all m > n. Since 
C"("-^)i/m {fi, 9) = Hn {(9) by Propositions we see that the polynomial Hn 
vanishes in a neig hborhood of 9o eR'^. Thus Hn ( x) is the zero polynomial 
and Hm {fi, 9) = for 9 G §^^^ and for all m > n. Hence we have proved b). 
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The implication b) — )■ c) follows from Kronecker's theorem. The implica- 
tion is c) — )■ a) is trivial. 

Clearly d) implies a), and it suffices to show that b) implies d). Now 
suppose that Hm (/^, ^) = for all m > n and that if„ (/i, 6) is not the zero 
function. Let Nn be the set of all 6 G S'^^^ such that Hn (/i, 6*) 7^ 0, so n is 
a normal index for each 6 G A^^. Clearly Nn is an open non-empty set. Let 
m > n he arbitrary. By Proposition 3.3 in [TBI P- 45] Tr^ (C; 6*) = 7r„ (C, ^) 
for all 9 E Nn and this implies vTn {(, 9) = Jl {(, 9) for all |C| > R and 9 G A^„. 
Let Pri (Cj 6') and (C? ^) as in the last section. Since the index n is normal 
for 9 G Nn it is clear that (P„ (C^ ^) ^ (C^ ^)) is an n-th Pade pair and we 
infer 

p„(c,^)/2(c,^) = g„(c,^). (22) 

By Theorem [7] there exist polynomials A{x) ,B (x) such that 

A m = C"Pn (C, 0) and CP (C^^) = C'Qn (C, ^) 

for all C and G S^^^. So we have for all y = p9 with p > R and 6' G A^„ that 

A{y)-fl,,^i{y) = \yfB{y) (23) 

where /ij-eai is defined in (jH]). Since /ireai {y), and obviously A (y) and B (y) , 
are real-analytic for all y G M'^, ||/| > P , the equality fl23|) . valid on an open 
subset of M°' \ Bji, holds for all y G M*^, \y\ > P as well. Moreover p^eai has an 
holomorphic continuation for all z E with L_ (z) > R where L_ is defined 
in on]). So we obtain 

A (z) pre^i (z) = {zl + ... + zl) B {z) for all z G C^, L_ {z) > P. (24) 

Suppose now that A contains an irreducible factor g which has a zero Zq G 
with L_ (2:9) > P. By continuity of L_ there exists a neighborhood f/ of Zq 
with L_ (2;) > P for all z G ?7. Equation (El]) shows that [/ H g^'^ {0} C 
P-^ {0} (recalling that zj + ... + zj for all z E with L_ (2) > P). 
It follows that g must divide P, see [23l p. 26]. Inductively, we can factor 
out each irreducible factor of A which has zero zq E with L_ (zq) > P. 
Finally we obtain polynomials Ai (z) and Pi (z) such that 

Ai (z) /ireai (^) = (-2? + ... + z^) Pi (;z) for all z eC^ with L_ (2) > P, 

and Ai (2;) 7^ for all L_ {zq) > P. The proof is accomplished. ■ 
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6 A cubature formula 



Let (T be a measure with finite moments on M and support in [—R, R] and 
consider the functional 

Tiu):=^[ «(C)?(C)c^C (25) 

where T^^ [t) = i?ie** for t G [0, 27r] for any Ri > R. For a polynomial 
u (C) = uq + ui( + ... + MmC™" "we have 

T{u) = Y,ui- fi = / u{x)da{x), (26) 

where /; := j^x^da {x) are the coefficients in the asymptotic expansion of a 
given in ([2]). 

We shall make use of the following classical fact (see e.g. [IB]): Let 
{Qni Pn) be the n-th Fade pair of the Markov transform a {() of a non- 
negative measure a with support in the interval [a, b]. If n is normal (so the 
Hankel determinant if^ (a) is not zero) then the zeros of Pn are 

real and simple and lie in the interval {a,b); moreover there exist positive 
coefficients such that the discrete measure 

cr„ = + ... + an5x„ 

is identical to a on the subspace of all polynomials p (x) of degree < 2?2 — 1 
and we have = Qn (xk) /P^ (a^fe) for k = 1, n. For any polynomial u (x) 
Cauchy's theorem yields 



o~ /" " p" dC = UkU (xfc) = I u{x) d(Tn. (27) 



Combining this with (1261) we obtain the following formula 

2n-l 

/ n (n ^^^dt = 

2m 



I «(oMv^c=Enz-/.. (28) 



valid for any polynomial u of degree < 2n — 1. 

Let V (M'^) be the set of all polynomials in d variables. 
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Definition 11 ^4 functional T : V (M'^) C is positive definite if 

T{u*u) > 

for all u (M*^) where u* is the polynomial obtained from u by conjugating 
the coefficients. 

Definition 12 A measure fi on with support in the closed ball is 
called Hankel-positive if the Hankel determinants are strictly positive, i. e. 

(l^, e)>0 for allneN,ee S'^'^ 

Obviously, an equivalent formulation for Hankel positivity is the require- 
ment that 

=0 1.. strictly positive definite 
for each 6 G E>'^~^. This means that for each 9 G S'^^^ and for all {xq, Xn) G 

W+\{xo,...,Xn)y^O 

n n 

h+i (^) ^i^i > 0- 

The following result is needed in the next theorem: 

Proposition 13 Let Pn{(,6) be defined in (f73]). If the Hankel determinant 
9 (—7- Hn (/X, 9) has no zeros then there exists Ri > such that 

((, ^) ^ for all 9 G QeC with |C| > Ri- (29) 

Proof. By assumption Hn {fJ',9) ^ for all 9 G S'^^^, so it follows that 
( I — )■ Pn (C, 0) defined in ( IT5i) is a polynomial of degree exactly n. Let us 
write 

Pn (C, 0) = Po (9) + p, (9) C + ... + Pn (0) C- 

Then p„ (9) ^ for all 9 G S'^^^ and Pn is continuous. A straighforward 
estimate now shows that there exists -Ri > such that P„ {(, 9) ^ for all 
Id > Ri and for all 9 G ■ 

The following is an analog of ( l26l) for the multivariate Markov transform, 
for the proof we refer to 
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Proposition 14 Let ^ be a measure on with support in and let Ri > 
R. Then for any u e P (M'^) 



M{u):= — — / / u{CO)iiiC,0)dCde = u{x)dfi{x). (30) 

The following result is an extension of the GauB quadrature formula to 
the multivariate setting. It can be seen as a solution of the truncated moment 
problem for the class of Hankel-positive measures. We refer to [1] and [9] for 
a description of the multivariate moment problem. 

Theorem 15 Let fi be a Hankel-positive measure with support in B^i. Let 
P„ (C) ^) (ind Qn (C) ^) be defined in l[T5\) and [T^ , and let Ri > R so large 
such that holds. Then the functional T„ : V (M*^) — C, defined by 

T^{u):=^—f f u{C0) Sj^^^dC de (31) 

for allu E V (M'^) , is positive definite and for each polynomial u (x) of degree 
<2n-l 

TAn) = ju(.)diAx). 

Moreover there exists a non-negative measure fin with support in an algebraic 
bounded set in M'^ such that 



Tn (u) = I U (x) (x) 

for any polynomial u. 



Proof. 1. Since (1291) holds for Ri > R the expression ( ]3T1) is well defined. 
Moreover {Qn {(, 9) , P„ (i^, 6)) is an ri-th Pade pair since Hn (/i, 6) ^ Q for all 
9 G S"^"^ by Hankel positivity. 

2. Suppose that m is a function of the form 

u m = uo (9) + ... + U2n-i (9) C'"-' (32) 
with continuous functions Uq, ...,U2n-i- By (!28|l we have 
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Integration over S'^ ^ gives 

^ 2n-l „ 

T„ (m) = — V / ui {6) fi {6) de. (34) 

3. Let u be a polynomial. Proposition [T^ shows that 

f u [x) dfi (x) = f [ u m l2 (C, 0) dCde. (35) 

If u has degree < 2n — 1 then for each 6 G E>'^^^ the function i— )■ m ((^^) is 
a polynomial of degree < 2n — 1, so m is of the form ( 1321) . Insert (jl]) in (1551) 
and integrate over F^^ to obtain 

2n-l „ 

uiie)fi{e)de. (36) 



/M (x) d/i (x) = — / 



Comparing with (15^ we conclude that T„ (w) is equal to f u (x) (i/i (x) 
for any polynomial of degree < 2n — 1. 

4. Let us discuss the question of positive definiteness of T„. Let R (x) be a 
real- valued polynomial. We have to show that T„ (R^) > 0. By the euchdean 
algorithm applied to the polynomials ( ^-^ R {(9) and ( Pn (C^) for each 
fixed 9, there exist a polynomial ( d{(,9) , and a polynomial C '"^ ^ (C^ ^) 
of degree < n, such that 

Rm = d{C,9)P^{C,9) + e{C,9). 

Write e (C, ^) = Cq (^) + ... + e„_i (9) ("-K Then 

(i? (C^))' = d' (C, ^) (P„ (C, ^)) ' + 2d (C, 9) e (C, ^^) Pn {C,9) + (C, ^^) . 



Multiply the last equation with Q„ {(, 9) jP^ (C, 9') and integrate with respect 
to C over r^^. Then 

Since C — ^ (C, 6') is a polynomial of degree < 2n — 1, fl33|) yields 

H0) = ^ [ (C, ^) I^T^c^C = E (^) (^) (^) • (38) 
"'rflj -r„ (C) 9) f^ i^Q 
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Integrate the last equation with respect to § ^ and use the definition of T„ 
in order to obtain 

1 r 

T„ [R^) =— ( V efc (^) Q {9) fk+i {e))dd. 
Js^-i ^ ^^^ 

Since is strictly positive definite we know that X]fc7=o ifi) ifi) fk+i [6) > 
for each 6 e in particular Tn{R^)>0. 

6. Let An be the polynomial defined in Theorem [7] such that An {(6) = 
C^^Pn (C^) • Note that An has real coefficients. It follows that for any poly- 
nomial u 

Tn (AnU) = / / C\ m Qn (C, 0) dC dO = 

since C, i — > C'^u {(9) Qn (C, 9) is a polynomial. The polynomial ( i — )■ P„ {(6) 
has only zeros in the interval {—R, R) , so it follows that P„ {p6) ^ for all 
p > R. Hence An {p9) ^ for all p > R, so the zero set of y i — > An (y) is 
contained in the ball B^. 

7. The existence of a representation measure pn follows from Theorem 
[16] below, cf. [19] (apphed to m = 2 and /i = An and /2 = —An). The proof 
is finished. ■ 

Theorem 16 (Schmudgen) Let & : V (W^) —^Cbea positive definite 
functional and let fi, fm be polynomials with real coefficients such that the 
set 

K := {xEW^ : fj{x) > for all j = l,2,...,m} 

is compact. Then there exists a measure p with support in K representing & 
if and only if 

&ifn-f,s-P*P)>0 (39) 
for all pairwise different ji, ...,js E {1, rn} and for all p eV (W^) . 

It is not difficult to see that a rotation invariant measure with non- 
algebraic support is Hankel-positive, cf . Theorem HI Now we give a different 
class of examples: 
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Proposition 17 Let Wq,Wi : [0, oo) — )■ M 6e bounded continuous functions 
with compact support such that \wi (r)| < Wq (r) for all r > 0, and Wq ^ 0. 
Assume that the measure fi has the density dfi := w (r, 'd) rdrdd where 

w (r, •&) = Wq (r) + Wi (r) cos-i? 

for all r > E [0, 27r] . T/ien /i Hankel-positive. 

Proof. Note that the assumption \wi (r)| < wq (r) for all r > assures 
that 

w (r, = Wq (r) + wi (r) cos ■& > (40) 

for all r > 0,^9 G [0, 27r] . Let us write 9 = e'"^ with e [0, 27r) . Note that 
Yo{e) = and Yfc,i(^) = -^cosA;^9 and 1^,2 (^) = ^^sinA;^?, k e No, 

provides an orthonormal basis of spherical harmonics. Then 

-| /-oo />27r /-oo 



|x| Iq (x) (jjU = / / r {r,d)rdrd'd = y2ii j r Wo{r)dr, 



and 



2 /'OO i'27v 

Ix]'^'^ Yi^i (x) d^ = I / r'^^ ■ r cost!} ■ w {r,{}) rdr d-d 

j^2s+2^^ (r) (ir, 

while all other distributed moments Cs,k,Tn are zero. By Theorem [T] we obtain 
the Markov transform: 

So /2s (e"') = r^'^+^wo (?^) and /2s+i (e*'') = cos-i? r'^^^'^wi (r) dr. 

Extend the function wq to an odd function Wq'^'^ on R \ {0}, so define 
Wq'^'^ (— r) := —Wq (r) for r > 0, and extend Wi to an even function w^^, so 

(— r) = t^i (r) for r > 0. Define a function : M — > M by 

(r) := r ■ [w"^'^ (r) + (r) cost?] . 

Note that (r) > for all r > by condition ( l40l) . Moreover (—''") = 
rwo (r) — rwi (r) cost? > for r > again by ( l40l) . Hence (r) > for all 
r G M. A straightforward calcuation shows that 

oo 



^ J -oo 
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for all / G No. This shows that {fi (e^'^))^^^^ is a positive definite sequence. 
If the sequence is not strictly positive definite then there exists a polynomial 
p (r) 7^ such that 




{p (r))^ Gi} (r) rdr = 0. 



Since G^ (r) is continuous on M \ {0} this implies that G^ (r) = for all 
r 7^ 0. Then = G^ (r) + G^ (— r) = 2rwo (r) for all r > 0, a contradiction 
to our assumption 7^ ■ 

In Theorem H] we have seen that the Markov transform /i {(, 6) of a ro- 
tation invariant measure has the property that ( 1 — > posseses an 
analytic continuation to the upper half plane. Next we show that the same 
is true for Hankel-positive measures. 

Theorem 18 Let ^ he a finite measure on with support in Bji. If is 
Hankel-positive then for each 9 G S'^""'^ the function C, 1 — > J1{C,6) posseses 
an analytic continuation to the upper half plane such that 

Imfi (C, ^) < for all ImC > 0, G S"'"^ 

Proof. Suppose that the sequence (// (6'));^q is strictly positive defi- 
nite. By the solution of the Hamburger moment problem (p. 65 in [18]) there 
exists a finite non-negative measure ae on M such that /i {(, 0) = j -^dag [t) . 
Hence ( 1 — ¥ fi (C, 0) extends to the upper half plane for C, and the condition 
Im Ji{(,6) < for all lm( > and 6 G S*^"^ follows from this integral 
representation. ■ 

Note that the measure ag in the last proof has the property that its 
support set is infinite since {fi{0))i^Qi is strictly positive definite. If we 
know that lmj2{(,9) < for all 9 G E>'^~^ and for all Im^ > then the 
function gg defined by gg {() '■= /^(C^^) is in the Nevanlinna class (see pQ) 
and the coefficients of the Laurent expansion are exactly the numbers /; (6) . 
Hence we can conclude that the sequence (/; (6*)); is positive definite for each 
9 G S"^^^. Note that Hankel positivity means that the sequence {fi (6*))^ is 
strictly positive definite for each 9 G S'^^^. 

Let us remark that in [15] we have introduced a different method for 
approximating a large class of signed measures, the so-called pseudo-positive 
measures, and we have provided a multivariate generalization of the Gauss- 
Jacobi quadrature formula. 

We thank the anonimous referee for the careful reading of the manuscript 
and for pointing to us a number of errors. 
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